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Abstract 

The formation of stripes in six-leg Hubbard ladders with cylindrical boundary conditions is investigated for two different hole 
dopings, where the amplitude of the hole density modulation is determined in the limits of vanishing DMRG truncation errors 
and infinitely long ladders. The results give strong evidence that stripes exist in the ground state of these systems for strong 
but not for weak Hubbard couplings. The doping dependence of these findings is analysed. 
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There is an ongoing controversial discussion about 
whether the ground state of interacting doped lattice mod- 
els in two dimensions like the t-J and the Hubbard model 
shows a charge modulation when subjected to particular, 
e.g., cylindrical boundary conditions. 

Recently, attention has turned to the two-dimensional 
Hubbard model on R x L-site ladders with local electron- 
electron repulsion U and electron hopping t, 
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where x = 1, . . . , R is the rung index and y = 1, . . . , L 
is the leg index. Cylindrical boundary conditions (closed in 
the rung y direction and open in the leg x direction) are 
assumed. For U = 0, this model describes a Fermi gas, ob- 
viously without stripes in the ground state. Moreover, us- 
ing renormalization group techniques, no instability toward 
stripe formation has been found in the weak-coupling limit 
U <C t [1]. In the strong-coupling limit U 3> t, the Hub- 
bard model can be mapped onto a t — J model with J = 
it 2 /U <C t, which exhibits stripes in the ground state, at 
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least on narrow ladders with J ~ 0.35t [2, and references 
therein]. An investigation of stripe formation in the Hubbard 
model at different couplings and fillings could thus signifi- 
cantly improve our understanding of these structures. 

In a recent density-matrix renormalization group 
(DMRG) calculation White and Scalapino [2] have shown 
that a narrow stripe appears in the ground state of a 7 x 6- 
site cluster for U > 6i. For weaker couplings, the hole and 
spin densities were interpreted as a broad stripe. However, 
no finite-size scaling has been performed, and the amplitude 
of the hole density modulation has not been investigated 
systematically as a function of DMRG truncation errors. 

In this work we exclusively consider 6-leg ladders (L — 
6) with R = 7r rungs for r = 1, . . . , 4. Since we are 
interested in the ground state of the hole-doped regime, we 
consider a system with Nx = 4r or N2 = 8r holes doped 
in the half-filled band, corresponding to RL — N\ = 38r 
or RL — N2 = 34r electrons. The average hole density is 
thus m = Ni/RL = 4/42 w 0.095 or n 2 = N 2 /RL = 
8/42 w 0.190, respectively. 

We employ a recently developed parallelised DMRG 
code [3], keeping up to m = 8000 density-matrix eigen- 
states per block for systems with up to R x L — 126 sites. 



We focus on the hole density h(x,y) 



>. 



where (...) represents the (DMRG) ground-state expecta- 
tion value. A stripe is related to a hole density modulation 
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Fig. 1. y-integrated hole density versus x for N\ — 4r (circles) and 
N 2 = 8r (squares) on a 21 X 6 ladder at U = lit. 

in the leg direction, 

L 

h(x) =^2h(x,y) . (2) 

y = l 

Stripe "signatures" derived from the staggered spin density 
are artifacts of the method and can be identified as such [4]. 

The study of stripe structures requires a spectral analysis 
of the hole density with respect to DMRG truncation errors 
and finite-size effects. The spectral transform is defined as 



H(k x , k y ) = J 1 x ^2 s'in(k x x)e' kyV h(x, y) (3) 

V x,y 

with k x — z x tt/(R + 1) for integers z x = 1, . . . ,R and 
k y — 2-KZy/L for integers — L/2 < z y < L/2. In the con- 
verged DMRG ground state we observe uniform behaviour 
of h(x, y) along the rungs. This implies that the spectral 
weight is concentrated at k y — 0. Stripes appear as hole 
concentrations which are translationally invariant along the 
rung (y) direction. At a hole doping of TVi = 4r, r stripes 
show up in a 7r x 6 ladder. When the doping is increased 
to N2 = 8r, the number of stripes doubles (see Fig. 1) and 
the structures become much less pronounced. 

The height and position of the maximum of the power 
spectrum [squared norm of Eq. (3)], i.e., the dominant har- 
monic can be extrapolated to the limit of vanishing DMRG 
truncation error (discarded weight W m )- This is possible 
because for small W m , expectation values of operators are 
polynomials in VW m . Indeed, we find a linear scaling of 
-ffmax = maxfc^. \H(k x , 0)| as \/Wm — > once the transi- 
tion to a striped state has occurred. For U = 12t, the ex- 
trapolated values of _ff ma x are finite. Thus we conclude that 
the hole density fluctuations found on finite ladders are not 
an artifact of DMRG truncation errors but a feature of the 
true ground state for U = 12t. For smaller values of U, 
the fluctuation amplitude decreases until U < 4t, and then 
stays close to zero if U is further reduced. 

Eq. (3) implies that when the limit 7? — > oo is taken, the 
amplitude of the dominant Fourier component in the hole 
density modulation (2) diverges linearly with \J~R. Conse- 
quently, results for infinite ladder length can be obtained by 
extrapolating H max /\fR as R^ 1 — > 0. Fig. 2 shows this 
limit for N — N\ (circles) and N — N2 (squares). 




Fig. 2. Amplitude of the dominant harmonic in the hole density modulation 
versus U, extrapolated to vanishing DMRG error and infinitely long ladders. 
Circles: N — 4r, squares: N — Sr. 

The following conclusions can be drawn from the data. 
For small U, stripe signatures observed in numerically de- 
termined ground states are artifacts of the method and vanish 
when proper extrapolation procedures are employed. Going 
from small to large U, there is a crossover from a homo- 
geneous to a striped state. With Ni = 4r, the transition 
occurs with a rather steep slope at U w At. Increasing the 
doping to N2 = 8r shifts the transition to larger U and 
makes it much smoother. Moreover, the number of stripes 
in the ground state is doubled. The existence of a similiar 
transition in real two-dimensional strongly correlated elec- 
tron systems would be of vital importance for the physics 
of layered high-T c cuprates. 
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